Slow conductance relcixation in insulating granular Al: evidence for screening effects 
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It is shown that the conductance relaxations observed in electrical field effect measurements 
on granular Al films are the sum of two contributions. One is sensitive to gate voltage changes 
and gives the already reported anomalous electrical field effect. The other one is independent of 
the gate voltage history and starts when the films are cooled down to low temperature. Their 
relative amplitude is strongly thickness dependent which demonstrates the existence of a finite 
screening length in our insulating films and allows its quantitative estimate (about lOnm at AK). 
This metallic-like screening should be taken into account in the electron glass models of disordered 
insulators. 

PACS numbers; 72.80.Ng, 61.20.Lc, 73.23.Hk. 



During the last 20 years, slow and glassy conductance 
relaxations were found in several disordered insulating 
systems [ll-0|. These relaxations have been revealed and 
often studied by means of field effect measurements. In 
such experiments, disordered insulating films are used as 
(weakly) conducting channels of MOSFET devices which 
allow us to measure their conductance response to gate 
voltage (Vg) changes. After a quench at e.g. AK, a 
slow and endless decrease of the conductance is found 
as long as Vg is kept constant, and any Vg change trig- 
gers a new conductance relaxation. Moreover, the system 
keeps some memory of its Vg history: Any stay under a 
fixed Vg value remains printed for some time in GiVg) 
sweeps as a conductance dip centered on this value. We 
have also shown recently in granular Al thin films that 
these Vg induced relaxations display aging, i.e. the dy- 
namics of the system depends on the time spent at low 
temperature, a characteristic property of glassy systems 

i- 

Several experimental findings suggest that these slow 
conductance relaxations could reflect the properties of 
the "electron glass" 0, 01j 8' gl^-^sy state theoretically 
predicted in the 1980s for disordered insulators Ac- 
cording to theoretical and numerical studies, a system of 
electrons with ill-screened interactions and disorder will 
need an infinite time to reach its equilibrium state at low 
temperature, the relaxation towards equilibrium being 
characterized by a conductance decrease of the system 
0. Up to now, the electron glass problem has received 
considerable theoretical developments but only few ex- 
perimental illustrations, which explain the interest in the 
electrical field effect results. Moreover, most of the ex- 
periments have focused on the electrical conductance and 
its V^-induccd relaxations but less is known concerning 
the dielectric properties and their possible time evolution 
at low T. 

Since even after a long stay under a fixed Vg the con- 
ductance decrease shows no sign of saturation, it is not 
possible to define the conductance relaxation relative to 
the equilibrium value, which is unknown. Instead, short 



excursions to Vg values never explored before are often 
used to define a Vg history-free reference conductance 
Gref 0-0| • This history- free conductance was sometimes 
called the "off-equilibrium" conductance since the system 
has never been allowed to equilibrate at these V^s. In the 
present Letter, we show that Gref is not constant in time 
after a cool down to 4:.2K in granular Al films thicker 
than lOnm. We also show that this feature demonstrates 
the existence of a screening length in our insulating gran- 
ular films and allows its quantitative estimate. 

Our granular Al films were prepared by e-gun evapo- 
ration of Al under a partial pressure of O2, as described 
elsewhere By changing the O2 pressure, we can tune 
the resistance of the films from metallic to insulating. 
For insulating samples, x-rays and TEM studies have re- 
vealed an assembly of crystalline Al grains with a typical 
size of a few nanometers. They are believed to be sep- 
arated by thin insulating AI2O3-X layers. MOSFET de- 
vices are made by deposition of granular Al films on top 
of heavily doped Si wafers (the gate) covered by a lOOnm 
thick thermally grown Si02 layer (the gate insulator). 

The film conductance G was measured by using a two 
terminal ac technique, employing a FEMTO current am- 
plifier DLPCA 200 and a lock in amplifier SR 7265. 
Source-drain voltage was such that G stays in the Ohmic 
regime. All the electronic equipments sensitive to room 
temperature drifts were placed in a thermalized chamber 
with a T stability better than O.IK . This last point was 
crucial in order to follow G variations with a precision of 
10"'^% during weeks of measurements. 

For 4:.2K measurements, the MOSFET devices are 
mounted in a box filled with He exchange gas and plunged 
into a 100^ liquid He dewar. The sample temperature 
follows the mK variations of the liquid He bath around 
4.2K and a carbon glass thermometer close by was used 
to correct these temperature variations. The time re- 
quired to cool down the sample from room temperature 
to 4:.2K is about lOmin. 

For Fig. [T]i, a MOSFET device with a 20nm thick 
granular Al channel was cooled down to 4.2K and main- 
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FIG. 1: (a) G{Vg) curves measured at different times t after 
a cool down to 4.27^. Vg = Vgsq — OV between Vg sweeps, 
(b) Corresponding G{t) curves for negative VgS in and out of 
the conductance dip region. The sample was 20nm thick with 
Rn = 550Mri at 4.2/C. See the text for the details. 



tained at this temperature under Vgeq = OV. Fast Vg 
sweeps from —15 to +15T^ (250s long) were taken every 
6000s after the cool down. All the G{Vg) curves display a 
conductance dip centered on Vgeq which reflects the mem- 
ory of the relaxation associated with the stay under Vgeq ■ 
Far enough from Vgeq (here for |V^| > 5V), G is roughly 
constant and this baseline reflects the "off-equilibrium" 
conductance Gref previously mentioned. Looking at the 
time evolution of G{Vg) curves, two features are salient. 
First, the amplitude of the dip increases as a function of 
time. This result is well known from previous studies 0- 
[^: the longer the stay under Vgeq, the more pronounced 
the dip is. Second, and this is the new feature we aim 
to discuss here, the baseline conductance also decreases 
as a function of time. As highlighted in Fig. [TJd, the 
G decrease is well described at any VgS by a hit depen- 
dence, more pronounced in the dip region because of the 
superposed baseline relaxation and dip growth. 

To be more quantitative, we can define a parameter 
whose physical meaning will become clear later, the slope 
ratio SR of the G relaxation slope in Fig. [TJd at Vgeq to 
the Gref relaxation slope (at —IQV for example). For 
the sample of Fig. [1] SR ~ 2.9. We have measured 4 



different 20nm thick films with Ro values from 10 Mil to 
lOGQ. A baseline relaxation was always observed, with 
SR. at AK between 2 and 3. For one film, the gate and 
the gate insulator were respectively Al and alumina and 
consistent results were obtained which excludes any role 
of a specific gate insulator material. Detailed investiga- 
tion of the vertical electrical homogeneity of the films has 
also shown that the baseline relaxation is not related to 
a specific layer close to the film-substrate interface . 

We checked that this baseline relaxation was indepen- 
dent of the Vg history by using the "two dip" protocol 
0. After some time under Vgeqi, Vgeq was changed to 
Vgeq2 for the rest of the experiment. We observed the 
formation of a new dip at Vgeq2 and the erasure of the 
old one at Vgeqi (as already known [1, 01), but the base- 
line relaxation in any Vg range far enough from Vgeqi and 
Vgeq2 continue as if no Vgeq change was imposed. 

One may think that the baseline relaxation is induced 
by the Vg sweeps themselves. Indeed, when measuring 
G{Vg) curves, a small dip starts to form at each measured 
point. The baseline thus reflects a short time relaxed 
G value (IDs being typical for our sweep parameters). 
We know from the "two dip" protocol results that a dip 
formed during 10s will be erased roughly in about 100s 
[3, H]. In our case, we wait for 6000s between two Vg 
sweeps, thus no memory of the previous sweep is expected 
out of the dip region. In one experiment, we checked that 
the baseline relaxation was unaffected by the suspension 
of Vg sweeps during 30ft-. 

An experimental artifact that may explain the baseline 
relaxation is the existence of a T drift after the cool down. 
Because of the rapid divergence of the films resistance at 
low T, the observed baseline variations correspond typi- 
cally to drifts of few niK (~ 2.5mK for the whole relax- 
ation observed at 4.2JC in Fig. [1]). We have performed 
different tests (thermometer stability, comparison with 
and without He exchange gas) which indicate that the T 
drifts are much smaller than the baseline relaxation am- 
plitude. The effect of the film thickness discussed below 
will definitely ruled out any interpretation in terms of a 
T drift. 

If the baseline relaxation is a property of the granular 
film itself, how can we explain it? A simple hypothesis 
is to state that the granular Al film, although electri- 
cally insulating, has a metallic-like screening length L^e 
that is smaller than its thickness (Th = 20nm in Fig. 
[1]). Then, only the layer of the film located at a distance 
smaller than Lgc from the gate insulator is sensitive to Vg 
changes. The conductance dip AGdip reflects the relax- 
ation of this layer, whereas the conductance baseline Gref 
reflects the relaxation of the rest of the film. The conduc- 
tance relaxation measured at Vgeq is then the sum of the 
two contributions: AG{Vgeq,t) = AGdip{t) + AGref{t)- 
Since the relaxations at any Vg are well described by a 
Int dependence, we can write AGrefit) = —Aref Ini and 
AGdip{t) = —Adiplnt. The slope ratio SR is by defini- 
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FIG. 2: G{t) measurements after a cool down at 4.2A' for: 
(a) a lOnm thick film, Ra = 6GQ; (b) a lOOnm thick film, 
Ra = 2M9.. The G{t) relaxations are plotted at W {Vgeq) 
and at — lOV^ (baseline relaxation Gref)- 



tion equal to {Adip + A^ef) I A^ef ■ If the film is homo- 
geneous, it is natural to suggest that Adip and Aref are 
respectively proportional to L^c (the Vg sensitive layer 
thickness) and Th — L^c (the Vg insensitive layer thick- 
ness). Then, SR = Th/iTh - L^c) as long as Th > L^c- 
The results for the 20nm films {SR bet"ween 2 and 3) 
give estimates for the screening length L^c bet"ween 10 
and 13nm. 

If our simple model is correct, then the baseline relax- 
ation relative to that measured at Vgeq must change "with 
the thickness of the films. This is illustrated in Fig. [2l 
"where the experiment of Fig. [1] "was reproduced for t"wo 
granular Al films lOnm and lOOnm thick. A clear thick- 
ness dependence is visible. Qualitatively, the baseline 
relaxation is almost absent for the lOnm thick film (Fig. 
[1^, SR = 20) "while it is almost equal to that measured at 
Vgeq for the lOOnm thick film (Fig. Hb, SR = 1.1). The 
agreement "with the above simple picture is also quan- 
titative. These SR values gives Lgc estimates of about 
lOnm, close to the results obtained from 20nm thick films 
(10 — 13nm). We stress that the SR variations observed 
for different thicknesses are much more important than 
those observed for a given thickness and different Rns. 

Our simple model also predicts a thickness depen- 
dence for the amplitude of the conductance dip AG/G = 
[Gref — G{Vgeq)]/Gref ■ AG/G data, mcasurcd for lOnm, 



20nm and lOOnm thick films are sho"wn in Fig. |3] as a 
function of their percolation critical resistance Rc [HI. 
In 2D films {Th less than the percolation radius Lq), 
Rc — Ro, but for thicker ones, Rc — {Th/Lo)Ra- Tak- 
ing Lq — 20nm [l2l, we estimate Rc as Rn for 10 and 
20nm films and 5Ra for lOOnm films. All the sam- 
ples "were measured after being kept under Veq = OV for 
t = 20h after a cool do"wn at 4:.2K. As can be seen and 
was already reported [l|,[3|, AG/G increases significantly 
with Rc- Besides, a clear reduction of AG/G is visible 
for a given Rc as Th increases, in qualitative agreement 
with the fact that the conductance dip originates from a 
layer of a fixed thickness. The agreement is also quanti- 
tative. As long as Th < Lgc, AG/G{t) does not depend 
on Th- Let us note its value C{t) (of course G{t) de- 
pends on Rc)- When Th > Lsc, one has AG/G{t) = 
Adiplnt/{Go - Arefhit) ~ Adiplnt/Go = C{t)L,e/Th 
{AG/G{t) < 1). Then, plotting {Th/L,e) x (AG/G) for 
all the samples should give one single curve G{t,Rc)- 
This is confirmed in the inset in Fig. [3] for which we have 
assumed a screening length of lOnm for all the films. 

Finally, our model provides a natural explanation for 
the increase of the AG/G data scatter with the thick- 
ness, as observed in Fig. [3l Small variations of the O2 
pressure or of the Al evaporation rate are always present 
during the deposition of granular Al films. Thus in thick 
samples, the measured Ru (and the estimated Rc) docs 
not necessarily refiect the resistance of the unscreened 
lOnm thick layer deposited on top of the gate insulator 
which determines the amplitude of the dip. The quality 
of the scaling in the inset in Fig. [3] is indeed an indirect 
test of the vertical homogeneity of the films, as discussed 
in details in Ref. (lo| . 

We now come to a discussion of our results. First, 
note that the thickness dependence of SR clearly demon- 
strates that the conductance relaxation exists in the 
"bulk"" of the film and is a property of the granular Al 
itself (and not of its interface with the gate or the surface 
oxide layer). The "bulk" relaxation after a cool down at 
A.2K is well described by a In t dependence without any 
sign of saturation over weeks of measurement. Such a In t 
dependence was also observed for the time increase of the 
conductance dip amplitude after a cool down [3, 0] (more 
complex laws are observed when Vgeq is fixed a time te 
after the cool down, see [1]). It is a natural relaxation 
law for a system having a 1/r distribution of relaxation 
times, as expected for an electron glass 13 1. 

In our granular Al thin films, the electronic transport 
at finite T is believed to result from electron tunneling 
between metallic Al grains. A lower bound for Lsc is 
then the typical grain size, which was found in recent 
TEM measurements to be 2 — 4nm. A second impor- 
tant microscopic length scale is the percolation radius 
Lq of the critical resistance network. Since in our films 
Rc < lOOGfi, the diffusion of charge carriers through the 
critical resistance network is expected to be faster than 
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FIG. 3: t^GjG values measured 20/i after a cool down at 4.27^ 
as a function of Rc (see the text for the details) . The films are 
lOnm (squares), 20nm (circles) and lOOnm (triangles) thick. 
The lines are guides for the eye. In the inset, (AG/G) x 
(Th/Lsc) versus Rc- Typical error bars on each point are 
about 10%. 

the time scale of the measurements (the intergrain ca- 
pacitance and the capacitance of the granular Al channel 
to the gate are respectively ~ 10~^^F and ~ lO^^^F). 
Thus, Lq should be an upper bound for L^c as indeed ob- 
served: Lq estimates from conductance fluctuations mea- 
surements are between 20nm and 40nm , a few times 
our Lsc value. 

We note that no baseline relaxation was reported on 
20nm thick indium oxide films which may indicate that 
Lsc{4:K) > 20nm in this system. Interestingly enough, 
Lq was found to be ~ 300nm in "crystalline" films jl4| . 
i.e. ~ 10 times our estimate for granular Al films, re- 
flecting a less dense percolating network in indium oxide 
consistent with a larger screening length. 

Our observations address the important question of 
screening in a disordered insulator, a subject which has 
been little explored both theoretically and experimen- 
tally. Electron glass models are generally developed 
in the limit of strongly localized electrons (Ref. 
[isl [lij are exceptions). But in real systems, there is 
a mobility of the charge carriers at finite T which will 
give rise to a metallic-like screening. Numerical stud- 
ies on disordered insulators have indeed found a transi- 
tion between a metallic screening at high T where most 
of the electrons are diffusive and a slow dielectric re- 
sponse at low T where most of the electrons remain lo- 
cated in finite size clusters (l7l - [l9j . Such a transition 
was observed around IK in capacitive measurements on 
a doped semiconductor According to the usual 

formula, L^^ = eeo/e'^dN/d^{EF), where dN/d^ is the 
thermodynamic DOS. In a glassy phase and for a finite 
time t, the system cannot explore all the configurations 



and the thermodynamic DOS has to be replaced by a 
pseudo-equilibrium DOS dN/dfj,{t) for which only relax- 
ations faster than t are allowed. The first theoretical 
attempts to do so 0, Hlj have found that the short time 
(~ Maxwell time) screening length di verg es as T goes to 
and decreases as a function of time [21|. 

In summary, we have described in insulating granular 
Al films the existence of a insensitive conductance 
relaxation. Its thickness dependence demonstrates the 
existence of a metallic screening length of about lOnm at 
AK. Our results provide a new way to study the screening 
length and its relaxation in disordered insulators where 
an anomalous electrical field effect has been found. They 
also point to the need for more theoretical studies: how 
are the classical electron glass models predictions affected 
by the existence of a metallic screening length ? 
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